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Abstract. 

We consider sequences of generalized Bell numbers B(n), n = 0,1,... which can 
be represented by Dobihski-type summation formulas, i.e. B(n) — ^(k) ' 

with P(k) a polynomial, D(k) a function of k and C = const. They include 
the standard Bell numbers (P{k) = k, D(k) = k\, C = e), their generalizations 
B r ,r{n), r = 2, 3, . . . appearing in the normal ordering of powers of boson monomi- 
als (P(fc) = D(k) = M, C = e), variants of "ordered" Bell numbers B ( f } (n) 
(P(k) = k, D(k) = { e y-) k , C = 1 + p, p=l,2...), etc. We demonstrate that for 

a, /3, 7, t positive integers (a, t ^ 0), [B(an 2 + [3n + 7)] * is the n-th moment of a pos- 
itive function on (0, 00) which is a weighted infinite sum of log-normal distributions. 
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In a recent investigation [1] we analysed sequences of integers which appear in the 
process of normal ordering of powers of monomials of boson creation a) and annihilation 
a operators, satisfying the commutation rule [a, a'] = 1. For r, s integers such that r > s, 
we define the generalized Stirling numbers of the second kind S r:S (n, k) as 



ns 



[(a^) r a s ] n = (a^ r -^J2 S rA^k)(^) k a k (1) 

k=s 

and the corresponding Bell numbers B r . yS {n) as 

ns 

B r , s (n)=J2 S r>,k). (2) 

k=s 

In [1] explicit and exact expressions for S rtS (n,k) and B r ^ s {n) were found. In a parallel 
study [2] it was demonstrated that B rs {n) can be considered as the n-th moment of 
a probability distribution on the positive half-axis. In addition, for every pair (r, s) 
the corresponding distribution can be explicitly written down. These distributions 
constitute the solutions of a family of Stieltjes moment problems, with B ryS (n) as 
moments. Of particular interest to us are the sequences with r — s, for which the 
following representation as an infinite series has been obtained: 

n— 1 

(3) 



(k 



k\ 



k=0 

l^ [k(k + l)...(k + r-l)] n 
eh (k + r-l)\ 



n > 0. (4) 

i /,' — / —II! 

k=0 

Eqs.(3) and (4) are generalizations of the celebrated Dobihski formula (r = 1) [3]: 

, oo , n 

5 i» = -Er n >0, (5) 

fc=0 

which expresses the conventional Bell numbers B lt i(n) as a rapidly convergent series. 
Its simplicity has inspired combinatorialists such as G.-C. Rota [4] and H.S. Wilf [5]. 
Eq.(5) has far-reaching implications in the theory of stochastic processes [6], [7], [8]. 

The probability distribution whose n-th moment is B r ^{n) is an infinite ensemble 
of weighted Dirac delta functions located at a specific set of integers (a so-called Dirac 
comb): 



o f n J 1 ^d(x-k(k + l )...(k + r-l I ^ n 

BrA*) = J o - j-E (jbT7=i)i p n ~°- (6) 

For r = 1 the discrete distribution of Eq.(6) is the weight function for the 
orthogonality relation for Charlier polynomials [9]. In contrast we emphasize that for 
r ^ s the B ryS {n) are moments of continuous distributions [2]. 

In this note we wish to point out an intimate relation between the formulas of Eqs. 
(3), (4), (5) and the log-normal distribution [10], [11]: 

. , 1 (InM-M) 2 

P^(x) = —^e 571—, x>0, <7,/i>0. (7) 



0.25 



e~ 1 8(x-1) 




x 

Figure 1. Weight function 1, l;x), see Eq.(13). 



First we quote the standard expression for its n-th moment: 

M n = r x n P^(x)dx = e n ^ +n ^) , n > 0, 
Jo 

which can be reparametrized for k > 1 as 



(8) 
(9) 



with 



t* = P Mk), (io) 

a = ^2a\n(k) > 0. (11) 

Given three integers a,/?, 7 (where a > 0), we wish to find a weight function 
Wi t i(a, (3, 7; x) > such that 

POO 

B 1A (an 2 + f3n + 7) = / ar n V^i,i(a, P, 7; x)da;. (12) 
Eqs.(5), (7) and (9) provide an immediate solution: 



(ln(a:)-/31n(fc)) s 



4a ln(fc) 



00 A; 7 exp 

6(x-l) + J2 . 

^ 2 2xk\\J 7raln(A;) 



(13) 



which is an infinite sum of weighted log-normal distributions supplemented by a single 
Dirac peak of weight e _1 located at x — 1. Thus it is a superposition of discrete and 
continuous distributions. Virtually the same approach can be adopted for the sequences 
B r>r {n), r > 1. In this case the k = 1 term in the numerator of Eq.(3) is larger than one 




Figure 2. Weight functions W r ,r(l> 1, 1; x) for r = 2, 3, 4. 



and so there will be no Dirac peak in the formula. Then the function W r ^ r (a, (3, 7; x) > 
defined by (a, (3, 7 integers, a, 7 > 0): 



/•oo 

.(an 2 + /3n + 7) = / x n W r>r (a, (3, 7; x)dx, 
Jo 



(14) 



is a purely continuous probability distribution given again by an infinite sum of weighted 
log-normal distributions: 



\ (k+r)\ 
I 00 [ k\ 



7-1 



cxp 



[ln(x)-/31n(^)] 2 
4aln[^ 



k=0 



2xk\\ net In 



(fc+r)! 



(15) 



The solutions of the moment problems of Eqs.(9), (12) and (14) are not unique. More 
general solutions may be obtained by the method of the inverse Mellin transform, see 
[12]. 

Several other types of combinatorial sequences have properties exemplified by 
Eqs.(12) and (14). We quote for example the so-called "ordered" Bell numbers B a (n) 
defined as [5]: 



B (n) = J2S(n, k)k\, 



(16) 



fc=i 



where S(n,k) are the Stirling numbers of the second kind, Si t i(n,k) in our notation. 
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These ordered Bell numbers satisfy the following Dobinski-type relation: 

i 00 un 

b ^) = - 2 Hy^ (17) 

fc=0 

from which a formula analogous to Eq.(13) readily follows. A more general identity of 
type (17) is [13]: 

1 00 , \ k n 

B<?\n) = — T ^P(^ T =J2S(n,k)k\p k , p=2,3,....(18) 

^ k = l / k = Q 

We will not discuss other types of sequences but rather observe that the relations 
of Eqs. (3), (4), (5), (17), (18) naturally imply that any power of these numbers also 
satisfies a Dobinski-type relation. As an example we give explicitly the simplest case of 
Eq.(5). For integer t > 0: 

™ r -KkJ%fc& <i9) 

with correspondingly more complicated formulas of a similar nature for powers of 
B rtr (n), B Q (n) and Bi P \n). For combinatorial applications of Eq.(19) see [14], [15] 
and [16]. 

We conclude that for any sequences of the type B(n) specified above 
[B(an 2 + (3n + 7)]* is always given as an ra-th moment of a positive function on (0, 00) 
expressible by sums of weighted log-normal distributions. We illustrate such a function 
for B h i(n) in Fig.(l). The application to B rr (n) for r = 2,3,4 is presented in Fig. (2). 
The area under every curve is equal to 1 on extrapolating to large x (not displayed). In 
both examples we have chosen a — f3 — 7 = 1. Observe the exceedingly slow decrease 
of these probabilities for x — > 00. This is confirmed by the fact that the moment 
sequences [B(an 2 + fin + 7)]' are extremely rapidly increasing. In the simplest case 
a = (3 = 7 = t = 1 we find B^n 2 + n + 1) = 1, 5, 877, 27644437, 474869816156751 for 
n = 0... 4. 

The circumstance that we can determine the positive solutions of the Stieltjes 
moment problem with both B(n) (discrete distribution) and [B(an 2 + (3n + 7)]' 
(continuous distribution) is a very specific consequence of the existence of Dobinski- 
type expansions. To our knowledge it has no equivalent in standard solutions of the 
moment problem. For instance, if the moments are n\ the solution e~ x does not give 
any indication as how one might obtain the solution for the moments equal to (n 2 )! . 

The strict positivity of W Tjr (a, ft, 7; x), for r = 1,2..., suggests their use in the 
construction of coherent states, which satisfy the resolution of identity property [17], 
[18], [19], [20]. This can be done by the substitution n\ — > B Ttr {an 2 + (3n + 7) in the 
definition of standard coherent states. More precisely for a complete and orthonormal 
set of wave functions \n) such that (n\n') = 5 n y and complex z we define the normalized 
coherent state as 

1 00 z n 

Mi> = 5 7K^TWTT) H - (20) 
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with the normalization 

which is a rapidly converging function of £ for < x < oo, x = \z\ 2 . Then, using 
the procedure of [18] we can demonstrate that the states of Eq.(20) along with Eq.(14) 
automatically satisfy the resolution of unity 

oo 

d 2 z\z;a,(3^}W r , r (a,(3,7; \z\ 2 )(z; a, (3, j\ = I = J2 l n )H, ( 22 ) 




n=0 



with 



W rA^(3,T,\z\)=n NMrM2) • (23) 

We are currently investigating the quantum-optical properties of states defined in 
Eq.(20). 

We close by quoting from [7] that, "the idea of representing the combinatorially 
defined numbers by an infinite sum or an integral, typically with a probabilistic 
interpretation, has proved to be a very fruitful one". In our particular case it has 
allowed us to reveal quite an unexpected relation between the Dobihski-type summation 
relations, which by themselves are reflections of boson statistics, and the log-normal 
distribution. 
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